Induced paramagnetic states by localized 7r-loops in grain boundaries 
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Recent experiments on high-temperature superconductors show paramagnetic behavior localized 
at grain boundaries (GB). This paramagnetism can be attributed to the presence unconventional 
d-wave induced 7r-junctions. By modeling the GB as an array of tt and conventional Josephson 
junction we determine the conditions of the occurrence of the paramagnetic behavior. 

PACS numbers: 74.72-h, 74.81. Fa, 75.20 -g 



The discovery of spontaneous currents in granular 
high-T c superconductors 0, 0, S S HI was a strong in- 
dication that a d-wave symmetry of the order param- 
eter is present in these materials. Indeed the d-wave 
scenario implies the possibility of existence of so-called 
7r-junctions, i.e., Josephson junction formed between su- 
perconductors with unconventional pairing which cause 
a 7r shift in the phase-current relation Q . A 7r-loop is an 
unconventional superconducting loop which contains an 
odd number of tt-j unctions. In zero field the ground state 
of a 7r-loop shows two energy degenerate magnetization 
states corresponding to two spontaneous current states, 
clockwise and counterclockwise. In non-zero magnetic 
field these spontaneous currents would act like orbital 
currents in paramagnetism Q. Therefore, if the sam- 
ple were field cooled, to permit the inner loops to feel 
the magnetic field, the response would be paramagnetic 
as indeed it was found in early works on Paramagnetic 
Meissner Effect (PME) in BSSCO by Braunisch et al. 

PME was observed also in different high-T c ceramic 
materials 0. 

However, the presence of PME in conventional low-T c 
samples shows that it cannot always be attributed 
to d-wave pairing []J| . Recently experiments and simula- 
tions were devised to test the relation between multiple- 
connectiveness and PME in conventional systems. A 
square array of LTC junctions was field cooled and shown 
to be paramagnetic over a large interval of magnetic field 
[Til H3L H^ . These papers also proposed a qualitative 
explanation for the effect based on the array multiple- 
connectiveness rather than the presence of tt-j unction. 
The effect of adding 7r-j unction in square arrays was an- 
alyzed in }l5| . 

The observation of spontaneous currents in YBCO 
biepitaxial 0°-90° tilt-tilt and twist-tilt grain boundaries 
(GB) junctions p indicates that paramagnetic effects 
due to d-wave pairing could be observed in GBs. In ad- 
dition, in Q spontaneous magnetic moments have been 
observed both in the High-T c films, where granularity or 
defects pin some vortices, and along the GB. Neverthe- 
less, the sample response in field cooling was diamagnetic. 
A recent experiment by E. IFichev et al. found that 
YBCO biepitaxial 45° asymmetric GB junctions in (nom- 
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FIG. 1: Mixed 71/ conventional one-dimensional Josephson 
junction array with localized ty- loops (half-gray). 



inally) zero field cooling show a paramagnetic response 
at low field. The origin of this paramagnetism could be 
debated. Is this simply due to presence of localized tt- 
loops that will act similarly to two-dimensional systems 
[lq or can be explained by means of paramagnetic quasi- 
particle currents due to existence of midgap states ? 
Here we want explore in detail the first alternative. 

In general a loop containing p Josephson junctions will 
have different magnetization states when a magnetic field 
is applied. If the junctions are identical the loop current 
I n is a solution of the following equation [lq : 
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where n = 0,l,..,p — 1 is the quantum number in the 
flux quantization expression. / is the frustation equal 
to the external flux normalized to flux quantum $0 and 
(3 is the SQUID parameter 2ttIoL/$>o with L the loop 
inductance and Iq the critical current of junctions in the 
loop. Varying n gives different families of independent 
solutions within a 2tt phase change k is an index 

which is equal to 1 if there are an odd number of ir- 
junctions in the loop and equal to zero otherwise. 

For any p the lowest energy solutions of Eq.(l) are 
diamagnetic for conventional loops and paramagnetic 
for 7r-loops H3. When (3 < 1 we have only one 
solution in the p = 1 loop which is diamagnetic in 
conventional loop and paramagnetic in 7r-loop without 
spontaneous currents 3]. But in multi-junctions loops 
(p > 1) we can have more states due to the presence 
of non-trivial solutions when changing quantum num- 
ber n. This implies that 7r-loops with, e.g., p = 2 
will show spontaneous currents also for low /?'s. Indeed 
for small (3 the solutions of Eq.(l) can be written as: 
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7± ~ sin(±7r/2-7r/)(l-cos(±7r/2-7r/)/3 L /2). So 
for / = we have two opposite spontaneous currents. For 
< / < 1/2 the solution 7 + is positive (paramagnetic) 
and 7_ is negative (diamagnetic). Moreover 7+ < 7 
giving a lower energy for the paramagnetic solution [2fl| . 

Small (3 7r-loops could be likely localized between GBs 
with different orientation along a junction 21] or where 
faceting cause an imperfect not completely flat GB pass- 
ing from a conventional junction to a tt-j unction or vicev- 
ersa. Recently also engineered "zigzag" arrays of mixed 
7r/conventional junctions have been realized and mea- 



sured (22, l22|. These can be described as an array of 
7r- loops separated by all conventional or all tt regions |24j . 

In the following we will describe the GB as an Id ar- 
ray of N + 1 Josephson junctions placed along it. The 
tt additional phase is supposed to vary along the array 
giving arise to tt and conventional sections separated by 
localized 7r-loops (see Fig.l) [13, Eij. We assume that 
system is not disordered. The magnetization dynamics 
of this TV-loops system can be described using the Dis- 
crete Sine-Gordon equation (DSG) |27j : 
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where cfj is the phase of the j-th junction in the GB, 

fj± — ^ e § j± is the frustation in the j^-th loop preced- 
ing (-) or following (+) the j-th junction; the index k(j) 
will be for conventional junctions and 1 for tt junctions. 
Times are normalized with respect to Josephson plasma 
frequency uoj and a is the normalized conductance. To 
include boundaries we set cpo = (fi, tpN+i = and 
fo — /iv+i = 0- We assume fj constant equal to / for 
1 < j < N. This implies that the magnetic field enters as 
boundary conditions on the two side loops of the array. 
The term ^4 is equal to the normalized magnetic field 
at boundary: 77 = ^ • \ L \jB ext (see Ref. 27]). Eq.(2) 
is analogous to that deduced in the continuous limit by 
E. Goldobin et al. in |24| in the context of analysis of 
"zigzag" arrays. We note that can be shown that Eq.(2) 
for TV equal one implies Eq.(l) for p equal two. 

In YBCO GB junctions the Josephson length Aj is 
smaller than GB physical dimension L, thus the normal- 
ized length I = L/Xj is larger than one The grain 
dimension along the GB Ax is usually smaller than L, 
being roughly of 1 jam for GB of Ref. 16] or also less in 
other circumstances |28j- GBs faceting is even smaller, 
ranging around 0.1 — 0.01 jam UtHl- A rough estimate 
of (3 can be made identifying (3 1 / 2 with the normalized 
length of grain ^ 0. From the data of Ref.s 0,0 is 
found A j ~ 5 fim which gives (3 c± 0.04. GBs faceting 
will give also a smaller (3. 

By integrating Eq.(2) we find the phases for all junc- 
tions. Initially the phases of conventional junctions are 
set to zero and the phases of tt-j unction to tt or — 7r, which 
are the stable equilibrium points of the single junction po- 
tential. This two possible choices correspond to two dif- 
ferent sign of the spontaneous current circulating around 
7r-loops. a = 1/yfPc was set to 0.25, which is within the 
interval proposed in Q. We do not use a field cooling 
process like in Q because initial conditions naturally set 



out diamagnetic or paramagnetic solution like in the sin- 
gle loop. In absence of bias current, the system naturally 
sets in a static equilibrium solution (ground state Q) 
after few plasma periods. Then the local magnetization 
is evaluated by: 
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where Acpj = ^+1 — (fj and the mean magnetization by: 
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where Aip = (£>n+i — <fi- In the absence of an exter- 
nal magnetic field the magnetization for a single local- 
ized 7r-loop in the array center (symmetric — tt junction 
[25^ have the shape reported in Fig. 2 topmost curves 
where the two spontaneous magnetization are shown for 
a N = 63 loop array with (3 = 0.04. The shape is very 
similar of "half-fluxon" obtained in the continuous ap- 
proach [2J due to relatively small (3. In Fig. 2 the effect 
of the magnetic field increase on the spontaneous mag- 
netizations is also shown. The magnetic field breaks the 
symmetry of two solutions: one is paramagnetic and the 
other diamagnetic. With the increase of the magnetic 
field the magnetization of the paramagnetic state is pro- 
gressively reduced due to the screening diamagnetic cur- 
rents that are generated at the boundary. The same cur- 
rents add to the magnetization of the diamagnetic state 
giving a larger diamagnetic magnetization. 

In Fig. 3 the mean magnetization for an array with a 
single 7r-loop is reported (circles). We note that magne- 
tization of paramagnetic state is zero at a threshold field 
77* c± 0.29. The linear decrease of mean is similar to that 
observed for (large (3) single loops 18]. For the parame- 
ters of Fig. 2 the physical threshold field is B* ~ 38 mG 
with the \l given in Ref. |j . 



3 



0.00- 






-0.01- 








-0.02- 




If 


1 


-0.03- 




1 


1 


-0.04- 






-0.05- 






-0.06- 


(a) f ; 




-0.07- 







10 20 30 40 50 60 
Section 



10 20 30 40 50 60 
Section 



0.05 
0.04 
0.03 
0.02 
0.01 
0.00 
-0.01 
0.02 



FIG. 2: Simulated magnetization of a N = 63 Josephson junc- 
tion array with a single 7r-loop in the middle with /3l = 0.04 
and a = 0.25: (a) diamagnetic solution with progressively in- 
creasing magnetic field r\ top to bottom 0, 0.1, 0.2, 0.3, 0.4, 
0.5 (b) paramagnetic solution with progressively increasing 
magnetic field [same values of (a)]. 



In Fig. 4a is reported the magnetization pattern in an 
array of N = 255 loops with 15 localized 7r-loops. Ac- 
cording to Ref. |24| flux quanta are sufficiently separated 
here to stay stable being the (minimum) length of conven- 
tional or 7r sections Ax/Xj ~ 4.64. The solution shows 
seven pairs positive-negative of half flux quanta plus an 
unpaired half flux quantum. In Fig. 4a the unpaired half 
flux quantum is positive, so solution is paramagnetic. An 
analogous diamagnetic solution exists when the unpaired 
half quantum is negative. Even 7r-loops configuration 
have zero spontaneous magnetization and are diamag- 
netic in small fields. Unpaired paramagnetic half flux 
quanta can be induced in the sample by a (moderate) 
field cooling process in small field, similar to 4]. The 
behavior of the mean magnetization is reported Fig. 3 . 
Both the spontaneous magnetization and threshold field 
are very small in this case. With the above data we find 
B* ~ 7.6 mG. In the same Fig. 4 is also reported the case 
in which 10 (Fig. 4b) and 12 (Fig. 4c) 7r-loops have ini- 
tial paramagnetic magnetization, which correspond to a 
stronger field cooling effect j20] . The corresponding mean 
magnetizations are again reported in Fig. 3. The mean 
magnetization for 12 paramagnetic 7r-loops becomes zero 
at 77* ~ 0.6 which corresponds to B* ~ 80 mG. 

For the sake of clarity and brevity, the results shown 
above have been obtained in absence of disorder. Disor- 
der has to be taken into account when we aim to describe 
high-T c materials and this will be subject of future inves- 
tigations. Here we just observe that disorder can locally 
change the penetration length altering the section length 
Ax/ X j and/or permitting larger screening currents in the 
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FIG. 3: Mean magnetization of both paramagnetic (upper 
curve) and diamagnetic (lower curve) solutions for Josephson 
junction mixed arrays. For all curves a — 0.25 and /3 — 0.04: 
o N = 63 with a single 7r-loop; N = 255 with 15 tt- loops and 
one odd paramagnetic half flux quantum; A N = 255 with 
15 7r-loops and 10 paramagnetic half flux quanta; □ N = 255 
with 15 7r-loops and 12 paramagnetic half flux quanta. 



sample. Small Ax/Xj implies that " currentless" (con- 
stant phase) states can occur 0, |2^j without sponta- 
neous currents. These facts, together with the small val- 
ues of the above threshold fields, imply that could be not 
surprising that also in moderate fields the state is dia- 
magnetic Q. Therefore, the presence (or the absence) 
of spontaneous currents would be no more strictly cor- 
related to paramagnetism. In ^| paramagnetism actu- 
ally appears without measurable spontaneous currents in 
Scanning SQUID microscope images. 

In conclusion localized 7r-loops in GBs can show both 
spontaneous magnetization and paramagnetic behavior. 
For samples large with respect to the penetration depth, 
implying a low (3 for each loop, paramagnetism exists in 
a relatively narrow region just near zero field. In absence 
of significant field cooling effects the energy difference be- 
tween diamagnetic and paramagnetic fundamental state 
solutions can be very small so observation of paramag- 
netism can be difficult or strictly depending on the par- 
ticular sample. Moreover in high-T c materials disorder 
can easily hinder the above picture. It is simpler to probe 
paramagnetic and diamagnetic states similar to that re- 
ported in Fig. 2 for engineered systems of 7r-loops as re- 
cently reported in Ref. |23] for two dimensional systems. 
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FIG. 4: Simulated magnetization of a N = 255 Josephson 
junction 15 7r-loop array with (3l = 0.04 and a = 0.25: (a) 
solution with one unpaired paramagnetic half flux quantum, 
top curve r\ = bottom curve r\ = 0.1; (b) solution with 10 
paramagnetic half flux quanta, top curve 77 = bottom curve 
77 = 0.5; (c) solution with 12 paramagnetic half flux quanta, 
top curve 77 = bottom curve 77 = 0.7. 
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